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Abstract 

Using stochastic quantization method [1], we derive gauge-invariant equa- 
tions, connecting multilocal vacuum correlators of nonperturbative field con- 
figurations, immersed into the quantum background. Three alternative meth- 
ods of stochastic regularization of these equations are suggested, and the cor- 
responding regularized propagators of a background field are obtained in the 
lowest order of perturbation theory. 



1 Introduction. 



Recently a new approach to investigation of the nonperturbative content 
of any field theory, based on the equations for vacuum correlators, derived via 
stochastic quantization method [1], was suggested [2]. It is closely connected 
to the Method of Vacuum Correlators [3], in which it is postulated, that the 
whole information about the QCD vacuum structure is maintained in the full 
set of irreducible vacuum averages (cumulants). 

Within this approach, applied to gauge theories, introducing correspon- 
ding generating functionals and using cumulant expansion [4,5], one obtains 
an infinite set of equations. These equations connect correlators, which con- 
tain various number of fields and stochastic Gaussian noise in gauge-invariant 
way. This approach is especially useful in the theories with nontrivial vacuum 
structure, e.g. in QCD, since the asymptotics of solutions of these equations 
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at Langevin time t tending to infinity yield the values of physical correlators 
without any assumptions about the model of the vacuum. In [2] the mini- 
mal closed set of such equations, corresponding to the Gaussian distribution 
of fields, was obtained for the case of gluodynamics and investigated in the 
lowest orders of perturbation theory both in standard, non-gauge-invariant, 
and gauge-invariant ways. 

The next related problem, one needs to solve in the framework of this 
approach, is the problem of separation of perturbative gluonic contributions 
in every term of cumulant expansion. To this end we split the total gluonic 
field into a background and a quantum fluctuation: 

A; = B;^gQ;, (1) 

where the principle of this division is unimportant [6]. In particular, the 
background field may be pure classical, and in this case we come to the 
problem of quantization of classical solutions, but, generally speaking, the 
fields form a quantum ensemble. 

Such a separation was used in [6-8] in order to develop perturbative theory 
in the confining QCD vacuum, which ensures the area law of an averaged 
Wilson loop with the value of the string tension, known phenomenologically, 
a ~ 0.2GeV^ [9]. In [7] this formalism was applied to the case of QCD at 
finite temperatures. 

It was shown in [6], that the confining background kills all the infra-red 
singularities, the lowest gluon and ghost corrections to the charge renorma- 
lization were calculated, and it was found out, that the usual logarithmic 
growth of as{R) in the empty space at large distances in the one-loop ap- 
proximation disappears in the presense of a background. Instead of that, one 
obtains, that as{R) is saturated at the scale of the inverse excitation mass of 
the transverse string vibration ~ 27r(7 ~ IGeV"^. 

In [8] unitary background gauges, where ghosts are either absent or non- 
propagating, were found, which may help to describe hybrid states in terms 
of physical polarizations of Q° only. 

However, in all the papers [6-8] the nonperturbative background fields 
were considered as given, and this input was parametrized, e.g. by the full 
set of cumulants. The main goal of this paper is to derive equations, starting 
from the Lagrangian, from which the correlators of background fields and of 
the quantum fluctuations may be obtained simultaneously. It seems na- 
tural to apply for this purpose stochastic quantization method [1], since the 
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background field formafism [10,6], developed within this type of quantization, 
possesses no-ghost property [11] as well as stochastic quantization of gluody- 
namics in the empty space [12]. The point is that all the individual stochastic 
diagrams, contributing to some gauge-invariant quantity may remain finite 
for t +00. The main idea [13] is to transform the gauge field in such a way, 
that in the Langevin equation for the new field the projector onto transverse 
degrees of freedom of gluonic field, standing in the action, is replaced by an 
invertible matrix. This transformation is required to leave unchanged all the 
gauge-invariant quantities, and, hence, should be a gauge transformation, 
but depending on t (because ^-independent gauge transformations leave the 
form of the Langevin equation invariant). 

It is known [12], that it is not necessary to add a gauge-fixing term into the 
Langevin equation, since the direct iteration in powers of coupling constant 
without introducing ghost fields leads to the same results as Faddeev-Popov 
perturbation theory, because the Langevin time takes the role of a gauge pa- 
rameter. For example, the coefficient at the projector onto longitudinal de- 
grees of freedom of the gauge field in the free propagator is lineary divergent 
at t tending to infinity, but, if one fixes calculates gauge-invariant quan- 
tities and then goes to the physical limit, t +oo, the divergent terms will 
cancel each other in the same manner as the terms, depending on the gauge 
parameter in the framework of the usual approach. That is why, in contrast to 
the Faddeev-Popov method of quantization, which reproduces correctly only 
small field fiuctuations (since for the case of strong fields Gribov ambiguities 
arise [14]), Langevin equation does not distort nonperturbative effects. 

Furthermore, it turns out, that the properly chosen ^-dependent gauge 
transformation modifies Langevin equation in such a way, that all the lin- 
eary divergent terms disappear completely, which is useful for calculations 
[13], and, in particular, suggests a method of quantization of non-holonomic 
systems [15]. This so-called stochastic or Zwanziger gauge fixing procedure 
was applied in [11] to quantization of gluodynamics in a background, and the 
/5-function in the one-loop approximation was computed. In what follows 
we shall exploit the Zwanziger term, introduced in the paper [11], to fix the 
gauge of a quantum fiuctuation, but, in contrast to [11], we shall not split 
Langevin equation into two parts in the sense of loop expansion, but use the 
total one in order to derive equations for correlators both of a background 
and of quantum fiuctuations. This work is performed in section 2. 

In section 3 we present three approaches to stochastic regularization of 
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the obtained equations. The first two covariant derivative approaches use 
the methods, suggested in the papers [16-19], while the alternative to them, 
but also Markovian, third approach is a new one. It is then shown, that all the 
three types of regularization leads to the properly regularized expressions for 
the propagator of a background field gluon in the lowest order of perturbation 
theory, when one neglects perturbative gluonic interactions. The main results 
of the paper and possible future developments are discussed in the Conclusion. 

2 Equations for correlators in bilocal approximation. 

In this section we present a general method of derivation of an infinite 
system of equations for correlators of background fields 5^, quantum fluc- 
tuations and stochastic noise fields rj^ and use it to obtain a minimal 
closed set of such equations, corresponding to the so-called bilocal approxi- 
mation, which follows from the assumption, that the grand ensemble of fields 
is Gaussian, so that all the cumulants, higher than quadratic, are put equal 
to zero. Lattice data suggest that this approximation has good accuracy in 
the confining regime of an averaged Wilson loop (for a discussion see the last 
reference in [3] ) . This hypothesis about the predominancy of bilocal correla- 
tions in the vacuum leads to the two alternative methods of investigation of 
higher correlators: the first one is based on the exact equations, where bilocal 
and higher correlators are considered on the same footing, while the other is 
the iterative one, where the values of bilocal correlators, obtained from the 
minimal system of equations, are then used to calculate threelocal correla- 
tors and so on. Moreover, in what follows we shall neglect all the quantities 
higher than of the second order in coupling constant, which means, that the 
equations to be obtained will not reproduce correctly those of Feynman dia- 
grams, which contain three- and four-perturbative-gluonic vertices. In order 
to extract explicitly the dependence on the coupling constant, we shall deal 
below with the fields = gA^, 6° = gB^ and = gQ^. 

A known important property of the background field method [10,6,8] is 
that it is possible to introduce the gauge fixing term for g^, which breaks down 
the gauge invariance of the partition function under quantum gauge trans- 
formations (which leave the background unchanged), but preserves gauge 
invariance under the so-called background gauge transformations 
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b,^u+ib, + id,)u, 

q, ^ U+q,U. (2) 

This background gauge condition leads to the unique choice of Zwanziger 
term, ensuring locaUty and renormahzabiUty of the theory [11], so that the 
Langevin equation takes the form 

. „ ^ (^W^W),. + g(D^^}D<Xr - ml, (3) 

where 



< 7]l{x, t)f]tix', t') >= 25^J^'5{x - x')6{t - t'), (4) 

and the sign of r/^ is changed. 

Due to (2), all the correlators, containing q*^ will be gauge-invariant, while 
for the background field one should use Schwinger gauge t) (x — xqJ^ = 
(where Xq is an arbitrary point), in which 6^ may be explicitly expressed 
through J^j^J: 



b^{x,t) = j dzya[z,x)T^^{z,t), 



where aiz.x) = ^^—Eoh q^^^ here and later in all the integrals of the type / 

the path of integration is a straight line. However, the final equations will be 
gauge-invariant in the same way, as it was discussed in [2]. 
Introducing the generating functional 



$/3 = Pexp ijdxy,< j dzya{z, x)J^l)'i}(z, t) + t) + (3 j dt'r],^(x, t') 



c 



(5) 

where C is some fixed closed contour and /3 is a c-number, one obtains, using 
Langevin equation (3): 
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c 

Applying to (6) the formula [5] 



(6) 



< e-^B >=< e'^>[<B> + Y.—,<- A"B » 



1=1 n'. 



(7) 



where A and B are two arbitrary operators, we have in bilocal approximation: 



d ( 
tr— < V^{y,xo,t)V^{u,xo,t) >= tr < V^{y,xo,t)lD^^\T^^J{u,XQ,t)^ 



1. d 
2 



^gSxf,Df^qp{u,xo,t))] >, 



(8) 



where 



^m(?/'^o,^) =/ dz^a{z,x)T^^^{z,XQ,t) + g\qp{y,XQ,t) + j dt'T]^{y,XQ,t,t') 



= ^(2^0, t)qi,{y, t)^{y, xo, t), 77;,(2/, xq, t, t') = $(xo, t)77;x(?/, t')^{y, xq, t), 



^{z, xq, t) = Pexp 
Noticing, that, due to (1), 



i j dz'^ha{z' ,t) 



one gets from (8) the first equation of bilocal approximation: 
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d \ 
+ < xo, t)Gpx^,{u, xq, t) >J , (9) 



where 



and we put all the terms with space-time derivatives to the right hand side. 
Here in order to disentangle the averages, containing covariant derivatives, 
we used the formulae: 

{d " 
-^M{u,xq)N + 1 1 dzaa{z,u)- 

■ {m{u, x^)NT^^}[z, Xo, t) - M{u, xq)T^^}[z, xo, t)N] 1 , (10) 



tr{DfD^^M{u, xo)N) = tr \ ^^q^M{u, xo)N + i(M{u, xo)NJ^^'^{ u,xo,t)- 



-M{u,xo)Tj)^^iu,xo,t)N^+i J dza[a{z,u)-^[M{u,xo)NJ^^^Xz.xo,t)- 



Xo 



-M{u,XQ)T^^}{z,XQ,t)N^ + ^a{z,u)(^M{u,XQ)NTl^}{z,XQ,t)- 

j + j dzaa{z,u) j dW(^a{w,u)- 

' Xo Xo 

'm{u, xo)J^jSiz, Xo, t)NJ^l^^{w, Xo, t)+M{u, xo)J^i^^{w, xo, t)NT^^}{z, xo, t)- 



-M{u, x^)T^^l {w, Xo, t)T^^}{z, Xo, t)N- 



-M{u, xo)NT^^} {z, Xo, t)J^i'^{w, Xo, t)jj, (11) 

where M(u, xq) is equal to J^j^^ (u, xq, t) or qa{u, xq, t), N = 
= N^^.,,^^{xi, ti, a:„, Xq, t) is, generally speaking, a product of some num- 
ber of Fa(3, Qa and ija, which are given in the points xi ^ u, ...,Xn ^ u at 
the moments ti,...,tn of Langevin time respectively, where all the parallel 
transporters between xq and each of these points are built of the field 6^ and 
given at the same moment t. 

Differentiating equation (6) twice by P, putting then P equal to 1, using 
the formulae (7), (10) and (11) and the definitions of three- and fourlocal 
path-ordered cumulants [4,2], one obtains two more equations of bilocal ap- 
proximation, where all the perturbative correlators, higher than of the second 
order in coupling constant are neglected: 

r d 

^^idi^ K(?/,a:o,t)/7/.(w,a;o,t,0 > - 
— fdv^f < V^{v, Xo, t)Va{w, Xo, t) > j < K(?/, Xo, t)r]^{u, xq, t, t') > + 



c c 

u 



+g J dz„a{z,u) j dw^{w,u)l< G pXfj,{u,xo,t)Tf^{w,Xo,t) > 

Xo Xo ^ 

< Tf}{z,xo,t)T]y{y,xo,t,t') > + 



+ < GpXi_iiu,xo,t)T]„{y,xo,t,t') >< T'i}{z,xo,t)T^pl{w,xo,t) > - 



- < Gpx^{u,xo,t)Tf}{z,xo,t) >< T]^{y,xo,t,t')Tfl{w,xo,t) > - 



- < GpXt,{u,xo,t)Tf^[w,xo,t) X r]y{y,xo,t,t')Tf}{z,xo,t) >j + 

u , 

+0^1 dzaa{z,u){< GpXp,{u,xo,t)Tf^{z,xo,t) >< r]y{y,xo,t,t')qx{u,xo,t) > + 



Xo 
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- < GpX^,{u,XQ,t)T]y{y,XQ,t,t') X qx{u,XQ,t)T^p}{z,XQ,t) > - 



- < Gpx^,{u,XQ,t)Tf}[z,XQ,t) X qx{u,XQ,t)r]y{y,XQ,t,t') >)} = 



^2 

= \ 9^ — < ^0, t)vAy, ^0, t, t') > + 



duxdux 


< -^v ^0, xo, t, t') > 



tr < 



J dt"-^ < riu{y, xo, t, t")r]^{u, xq, t, t') > + < r]y[y, xq, t, t)r]^[u, xq, t, t' 



dt 



- < r]„{y, xo, t, t')r}p{u, xq, t,t) > - 

\jdv^j dWa < V^{V, Xq, t)Vaiw, Xq, t) > 



c c 



t / u 



• j dt" < T]^{y,xo,t,t")T]ij,{u,XQ,t,t') > -\- j> dza j dt"ij dw(^a{w,u)- 

C a^o 

•(< •^A2(^'^o,^)^a(^,a;o,t,0 ><r]u[y,XQ,t,t')Tfl[w,XQ,t) > - 

- < -^Am ^0, t)Tf^{w, xo, t) X r]aiz, xo, t, t")r]y{y, xq, t, t') + 

+§{< Gpx,jiu,xo,t)r]a{z,xo,t,f) >< r]i,{y,XQ,t,t')T^p^[u,XQ,t) > + 



9 fdzaj dt"i j dW(^(a{w,u)-^(< Gpx^{u,XQ,t)Tf^{w,XQ,t) > 



■ < r]y{y,XQ,t,t')T''xQ{w,XQ,t) >j + ^a{w,u)(^ GpXp{u,XQ,t)- 

■^pciw,Xo,t) >< T]^{z,Xo,t,t")T]^{y,Xo,t,t') > - 

< GpXp[u,XQ,t)T]a[z,XQ,t,t") ><T]y[y,XQ,t,t')Tfl[w,XQ,t) + 



^^(^ < Gpx,,{u,XQ,t)r]a{z,XQ,t,t") >^ < r]y{y,XQ,t,t')qx{u,XQ,t) > - 



-g < xo, t)r]a[z, xo, t, O > ^ < 77i.(y, xo, t, t')GpXp{u, xo, t) >j . (13) 

As was discussed in [2], in the physical hmit, t +oo, in the confi- 
ning regime of an averaged Wilson loop, the dependence on the point xq is 
negligible, since the difference between each of the cumulants and its gauge- 
invariant analog is of the order of < 0.04 [2], where Tg is the correlation 
length of the vacuum, at which a cumulant vanishes and is a space width 
of a Wilson loop, and we obtain gauge-invariant equations for correlators of 
the fields T^^} and q^^ 7]^, immersed into a background as insertions in the 
parallel transporters. 

Note, that among these equations only the first one, equation (9), is linear, 
while equations (12) and (13) produce complicated hierarchy of perturbative 
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correlators up to the second order of perturbation theory. Hence, the pro- 
cedure of solution of equations (9), (12) and (13) is the following: first one 
needs to put all the perturbative fluctuations equal to zero and to solve the 
so-reduced equations for correlators of a background and of the Gaussian 
noise fields (which are just the equations (17), (19) and (20) from the paper 
[2]). After that one should to include perturbative interactions, expanding 
the correlators, containing q^j^ up to the order of and using the obtained 
values of pure background and noise fields' correlators. 



3 Stochastic regular izat ion and perturbative expan- 
sion of equations (9), (12) and (13). 

In this section we present three methods of stochastic regularization of 
equations (9), (12) and (13) and use them to derive regularized propagators 
of a background in the lowest order of perturbation theory in the absence of 
perturbative corrections. The first two of them are based on covariant deriva- 
tive regularization schemes, which were suggested in the papers [16-18] and 
used in [19] to calculate /3-function in QCD in the one-loop approximation. 

First is the so-called power-law regularization scheme 



dyR'\x,y,t)T]l{y,t)., 



where 



1 



-\ ah 



{x,y,t), 



n=l,2, 



A is an ultraviolet cutoff, /Y°{x,y,t) = J dz{D^Y'^{x, z,t){D^Y°{z,y,t) is the 
covariant Laplacian with {D^Y^{x,y,t) = D^^^"-^(x,t)S{x — y). One of the 
main results of the paper [17] is the proof, that any Yang-Mills theory in d 
dimensions is regularized to all the orders, when we choose n > [^{d -\- 1)], 
where [x] is the largest integer less than or equal to x. 

The generalization of equations (9), (12) and (13) after applying such a 
regularization is obvious. For example, equation (9), written in details, takes 
the form: 

2^^qI]J dz\[z,y) j dxpa{x,u) < T^^J{z,Xo,t)T^^J{x,Xo,t) > + 
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/ 

+p J dzxa{z,y)l<J^^J(z,xo,t)q,,(u,xo,t) 



> + 



-\- J dt' J dw < Tx}{z^XQ^t)^i^{u^w^XQ^t^t') >\+g J dxpa{x,u)i<q„{y,xo,t)- 

^ Xo ^ 



•TfJ{x,xo,t) > + j dt' j dw < ^,y{y,w,xo,t,t')T^^J{x,xo,t) 



> + 



+9^ < (lv{y,XQ,t)q^{u,XQ,t) > j dt' j dw- 



< q^{y, Xq, t)^f,{u, w, Xo, t, t') > + < ^^{y, w, Xq, t, t')q^{u, Xq, t) > + 



+ j dt" Jdv < $,j,{y, w, xq, t, t')^^{u, v, xq, t, t") > ) [ = 



= tr 



d 



dux 



d d 
+^/-— < T^}{z,XQ,t)q^{u,XQ,t) > -g^— < T^}{z,XQ,t)qx{u,XQ,t) > + 



dux 



du, 



+fi'^A/x^ < J^i^Jiz,XQ,t)qp{u,xo,t) >j+g< qi,{y,X(i,t)Tfl{u,XQ,t) > + 
t 

+g j dt' j dw < ^u{y,w,XQ,t,t')Txl{u,XQ,t) > + 



d 



d 



+5' 1 ^ < (li^iy^ xq, t)q^{u, Xo, t) > < q^{y, xq, t)qx{u, xq, t) > + 



Mux 



d 



+^A/x^ < qv{y,XQ,t)qp{u,XQ,t) > + 
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+ j dt' j dw(y-^^ < ^^{y,w,xo,t,t')q^{u,xo,t) > - 
d 



d 



< iy{y, w, xq, t, t')qx{u, xo, t) > + 



+(^A/x^ < ^y{y, W, Xq, t, t')qp{u, Xq, t) > 



where ^Ju^w^xo^t^t 



p 

/\ 



^{xo,u,t) 



ab 



M, t')rli{w^ f)t'^^{w, Xo, t). 



The equations (12) and (13) changes correspondingly. 

In the lowest order of perturbation theory, when one neglects perturbative 
fluctuations, and equations (9), (12) and (13) reduce to the equations for 
background fields only, one may check in the same way, as it was done in [2] 
for the unregularized case, that for n = 2{d = 4) the regularized propagator, 
obtained from these equations, has the form: 



+ 2kHL^X (14) 

where T^j, = 6^1^ — ^^^^Ljj^jj = are the transverse and the longitudinal 
projectors respectively. In (14) one may recognize the expression (2.19) from 
the paper [17], where the "gauge-fixing parameter" a = 2kH, which seems 
to be natural in the sense, stated in the Introduction. 

Second method of regularization exploits the so-called heat-kernel re- 
gularization scheme [18], where the regulator has the form R"'^{x,y,t) = 
= (e^y^{x, y, t), that is believed to be technically superior for nonperturba- 
tive analysis. In analogous way, one gets from the equations (9) and (12) in 
the lowest order, when the perturbative gluons' contributions are neglected: 

(15) 

which coincides with the formula (9) from the paper [18] oX a = 2kH. 



13 



Finally, let us present a new one, also Markovian and preserving gauge 
invariance, type of regularization |^. Its basic idea is to smear space-time 
delta-function in (4) so, that properly modified Langevin equation remains 
gauge-invariant. To this end we introduce the new noise fields 



i;{x,t) ^^-j-tjdy e ^v;{y,t), (16) 



47^2 



SO that 



< C,{x,mt{y,t') >= ^^"''^(^ - t')e-'^, (17) 



hm < e,{x,t)C{y,t') >= 26,J-'6{x - y)6{t - t') 

A|— ^ + CXD 

and modify the Langevin equation (3) in gauge-invariant way 



/ A 2\2 2 2 

- / dy e~^^^^x,y,t)r]^{y,t)^y,x,t), (18) 

keeping in mind, that at | A |^ +oo the integral in the right hand side is 
saturated at | 7/ — x |<C j^. 

Then, due to (16) and (17), in the lowest order of perturbation theory 
in the absence of perturbative gluons, it follows from the equations (9) and 
(12), regularized according to (18), correspondingly 



^This method was suggested by Professor Yu.A.Simonov. 
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< Bl(u,t)Bl(y,t) > +^ / dt' j dve-^^ < Bl{u,t)r,t{v,f) >j, 

(19) 







■ J dve-^^ < Bl{u,t)Ti[v,t') >= - t')e-^^. (20) 

Looking for < B'j^(y,t)r]^^{u,t') > in the form 6"'^d^iy{z,T), where z = u — 
-y,T =\ t -t' \,d^^{z,T) = d^^{z,T), d^^{-z,T) = d^^{z,T), one obtains 
from (20): 

d,u{k, t) = -^^(T)e-A-(T^.e-^ ^ + V), (21) 

where J^j.(a;,r) = /d?/e ^^o?^j.(?/, r). 

Looking for < i?°(^, t)B^^{u, t) > in the form 5^^h^y{z^ t), where hiy^{z, t) = 
= hi^iy{z,t),hn^{—z,t) = h^j,y{z^t) and using (4) and (20), one gets from the 
equation (19): 

and, hence, due to (21), the regularized propagator of a background field has 
the form 



(22) 

Therefore, we see, that this method of regularization is similar to the heat- 
kernel method, but is simpler, since there do not arise higher derivatives in 
the regulator. 
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4 Conclusion 



In this paper we applied stochastic quantization [1] to develop a method 
of derivation of an infinite system of exact equations for gauge-invariant 
correlators in gluodynamics, where all the perturbative contributions are ex- 
tracted explicitly in the form of insertions into background parallel trans- 
porters. Therefore, the obtained equations allow one to derive pure back- 
ground correlators and the correlators, containing perturbative corrections, 
simultaneously, using for quantization the same stochastic noise fields. 

After that we obtained the minimal set of equations of bilocal approxima- 
tion (corresponding to the Gaussian distribution of fields), where we threw 
away all the perturbative interactions higher than of the second order, and 
suggested for it three methods of stochastic regularization, all of which pre- 
serve gauge invariance of the obtained equations. The first two of them, based 
on the so-called covariant derivative regularization schemes, lead, in the low- 
est order of perturbative theory in the absence of perturbative gluons, to the 
known values of the regularized propagator of a background field, while the 
third method is a new one. It yields the results, similar to the heat-kernel 
method, but is simpler than the latter, since in the framework of this method 
higher derivatives in the regulator do not exist. 

The application of the suggested approach to treating the large- A/^ regime 
of QCD as well as the new equation for the Master field and its connection 
with the Bootstrap equation will be a topic of a separate publication. Possible 
types of solutions of the derived equations will be presented elsewhere. 

The results, presented in this paper, were partially reported at the In- 
ternational Workshop " Nonperturbative Approaches to QCD", Trento, Italy, 
July 10-29, 1995. 
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